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Abstract 

We present a measurement of the time-dependent CP violation parameters in B° — > tt + tt~ 
decays. The results are obtained from the final data sample containing 772 x 10 6 BB pairs collected 
at the T(45) resonance with the Belle detector at the KEKB asymmetric-energy e + e~ collider. 
We obtain the CP violation parameters 

A C p(B° tt+O = +0.33 ± 0.06 (stat) ± 0.03 (syst), 
S C p(B° -> 7T+7T-) = -0.64 ± 0.08 (stat) ± 0.03 (syst), 

where Acp and Sep represent direct and mixing-induced CP violation, respectively. From an 
isospin analysis including results from other Belle measurements, we rule out 23.8° < (p2 < 66.8° 
at the la level. 
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FIG. 1: Leading-order diagrams contributing to B° — > tt + tt" decays, (a) shows the dominant 
lst-order (tree) while (b) shows the 2nd-order loop (penguin) diagram. In the penguin diagram, 
the subscript x in V x b refers to the flavor of the intermediate-state quark (x = u, c, t). 



INTRODUCTION 



Violation of the combined charge-parity symmetry (CP violation) in the standard model 
(SM) arises from a single irreducible complex phase in the Cabibbo-Kobayashi-Maskawa 
(CKM) quark- mixing matrix [j], 0|. The main objective of the Belle experiment at KEK, 
Japan and the BaBar experiment at SLAC, USA, is to over-constrain the unitarity triangle 
for B u ^ decays. This allows a test of the CKM mechanism for CP violation as well as the 
search for new physics (NP) effects beyond the SM. Mixing-induced CP violation in the B 
sector has been clearly observed by the Belle j3[ and BaBar 0] Collaborations in the b —> ccs 
induced decay B° — > J/ipK^, while many other modes provide additional information on CP 
violating parameters. This paper describes the measurement of CP violation parameters of 
B° — » 7t + tt~ decays, shown in Fig. [TJ 

Decays that proceed dominantly through the b — > uud transition are sensitive to the 
interior angle of the unitarity triangle 2 (a) = arg(— V t dV t * b ) / (V u dV* b ) . The Belle, BaBar 
and LHCb experiments have reported time-dependent CP asymmetries in these modes that 
include decays such as B° tt+tt" (pnf @, 8, P + P~ 0, EH and afn^ [H, [l3 . 

The decay of the T(4S) can produce a B°B° pair that must oscillate coherently, of which 
one (-Br 6C ) may be reconstructed in the 7r + n~ final state. This final state does not determine 
whether the B^ ec decayed as a B° or as a B°. The other B meson in the event {B^ &g ), 
however, can be reconstructed in a final state that determines its 6-flavor, and therefore the 
flavor of the B^ cc is opposite to the B^ ag flavor at the time of its decay. The proper time 
interval between -Br 6C and B^, which decay at time t^ ec and txag, respectively, is defined 
as At = tR ec — txag- For the case of coherent B°B° pairs, the time-dependent decay rate for 
a two-body decay when B^ possesses flavor q (B°: q = +1; B°: q = —1), is given by 



-|At|/r s0 



Acp cos ArridAt + Sep sin Am^At 



(1) 



Here, r^o is the B° lifetime and Arrid is the mass difference between the two mass eigenstates 
of the neutral B meson. This assumes CPT invariance, no CP violation in the mixing, 
and that the difference in decay rates between the two mass eigenstates is negligible. The 
parameter Acp measures the direct CP violation, while Sep measures the amount of mixing- 
induced CP violation. 
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FIG. 2: Complex isospin triangles from which A^>2 can be determined. 



A feature common to the B° — > tt + tt~ final state and the other b —> uud modes men- 
tioned earlier is that an effective angle <p^, is determined rather than $2 itself due to the 
possible presence of additional penguin contributions that may interfere with the dominant 
tree contribution. Thus, sensitivity to </> 2 comes from this relation between the measured 
parameters Sep = yl — A? CP sm 2<^>§ H . In the limit that only the dominant tree amplitude 
contributes, no flavor- dependent direct CP violation is expected and (j)^ = <p2- 

Despite penguin contamination, it is still possible to determine A02 = 0| ff — 02 in B° — > 
tt + tt~ with an SU(2) isospin analysis [lij by considering the set of B — > tttt decays into the 
three possible charge states for the pions. Here, the two pions in B + — > tt + tt° decays must 
have a total isospin of / = 1 or / = 2, since J 3 = 1. For the penguin contributions, only 
/ = or / = 1 is possible since the gluon is an isospin singlet carrying 1 = 0. However, 
I = 1 is forbidden by Bose-Einstein statistics, thus strong loop decays cannot contribute 
and hence B + — > 7r + 7r° decays only through the tree diagram in the limit of neglecting 
electroweak penguins. 

The complex B — > tttt decay amplitudes obey the relations 

A +0 = ^|^+- + A)o, ^-o = ^ A +~ + ^oo, (2) 

where the subscripts refer to the combination of the pion charges. The decay amplitudes 
can be represented as triangles shown in Fig. [2j As B + — > tt + tt° is a pure tree mode, 
these triangles share the same base, A + q = A^q, and A<p2 can be determined from the 
difference between the two triangles. These triangles and 02 can be fully determined from 
the branching fractions, B(B° tt + tt~), B{B° tt°tt°), B(B + tt + tt°) and the CP 
violation parameters, Acp(B° — > tt + tt~), Scp(B° — > tt + tt~) and Acp(B° — > tt°tt°). This 
method has an 8-fold discrete ambiguity in the determination of <p2, which arises from the 
4 triangle orientations about A +0 and the two solutions of (ftf^ in the measurement of Sep- 

The Belle, BaBar and LHCb Collaborations have reported measurements on all param- 
eters reported here in the past, which are collected in Table [H The results of the previous 
Belle measurement was based on a 535 million BB pair sample and are superseded by the 
analysis presented here. 

In Sec. (Til we briefly describe the data set and Belle detector. We explain the selection 
criteria used to obtain signal candidates and suppress backgrounds in Sec. 1111} followed by 
the fit method used to extract the signal component in Sec. [TV] In Sec. |VJ the results of the 
fit are presented along with a discussion of the systematic uncertainties in Sec. ED Finally, 
our conclusions are given in Sec. IVHl 
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TABLE I: Summary of published CP violation parameters obtained by the Belle j5|], BaBar [g] 
and LHCb Collaborations 0]. For all parameters, the first uncertainty is statistical and the sec- 
ond is systematic. The Belle value for Acp is marginally consistent with the BaBar and LHCb 
measurements. 



Parameter 


Belle Value 


BaBar Value 


LHCb Value 


A C P(B° -»• TT+1T-) 


+0.55 ±0.08 ±0.05 


±0.25 ±0.08 ±0.02 


+0.11 ±0.21 ±0.03 


S C p{B° -> 7T+7T-) 


-0.61 ±0.10 ±0.04 


-0.68 ±0.10 ±0.03 


-0.56 ±0.17 ±0.03 



II. DATA SET AND BELLE DETECTOR 

This measurement of the CP violation parameters in B° —> tt + 7t~ decays is based on 
the final data sample containing 772 x 10 6 BB pairs collected with the Belle detector at 
the KEKB asymmetric-energy e + e~ (3.5 on 8 GeV) collider (H)]. At the T(4S) resonance 
(y/s = 10.58 GeV), the Lorentz boost of the produced BB pairs was (3 / -f = 0.425 along the 
z direction, which is opposite the positron beam direction. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprising of CsI(Tl) crystals (ECL) located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [l6|. Two inner detector configurations 
were used. A 2.0 cm radius beampipe and a 3-layer silicon vertex detector (SVD1) was used 
for the first sample of 152 x 10 6 BB pairs, while a 1.5 cm radius beampipe, a 4-layer silicon 
detector (SVD2) and a small-cell inner drift chamber were used to record the remaining 
620 x 10 6 BB pairs We use a GEANT-based Monte Carlo (MC) simulation to model 



the response of the detector and determine its acceptance [18 



III. EVENT SELECTION 

The decay channel B° — y 7r + 7r~ is reconstructed from two oppositely charged tracks. 
Charged tracks are identified using a loose requirement of distance of closest approach to 
the interaction point (IP) along the beam direction, \dz\ < 4.0 cm, and in the transverse 
direction, dr < 0.4 cm. With information obtained from the CDC, ACC and TOF, particle 
identification (PID) is determined with the likelihood ratio Ci/j = £«/(£« + £j). Here, £j 
(Cj) is the likelihood that the particle is of type i (J). To suppress background from particle 
misidentification, vetoes are applied on particles consistent with the electron hypothesis. 
The PID ratios of the two charged tracks are used in the fit model to discriminate 

between the 3 possible peaking channels that may occupy the final state: Bz — > K~7i + 
and B° — y K + K~ . The additional SVD requirements of two z hits and one r — hit 



are imposed on the charged tracks so that a good quality vertex of the reconstructed B 
candidate can be determined. 

Reconstructed B candidates are identified with two nearly uncorrleated kinematic vari- 
ables: the beam-energy-constrained mass M^ c = \f (-^beam) 2 — (Pb MS ) 2 , and the energy dif- 
ference AE = Eb MS — E™^, where E™^ is the beam energy and E^ MS (p§ MS ) i s tne energy 
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(momentum) of the B meson, all evaluated in the e + e~ center-of-mass system (CMS). The B 
candidates that satisfy M bc > 5.24 GeV/c 2 and -0.20 GeV < AE < 0.15 GeV are selected 
for further analysis. 

The dominant background in the reconstruction of B^ cc is from continuum (e + e~ — > qq 
where q = u,d, s, c) events. Since their topology tends to be jet-like, in contrast to the spher- 
ical BB decay, continuum events can be distinguished from BB events using event-shape 
variables, which we combine into a Fisher discriminant J^^i/qq 20[. The BB training sample 



is taken from signal MC, while the qq training sample comes from data taken below the 
T(4S) resonance. The Fisher discriminant is then constructed from the variables described 
in Ref. [l2[ • The variable providing the strongest discrimination against continuum is the an- 
gle between the B^ cc thrust direction (TB) and the thrust of the tag side (TO) | cos^tb^oI- 
The thrust is defined as the vector that maximizes the sum of the longitudinal momenta of 
the particles. For a BB event, the pair is essentially at rest in the CMS, so the thrust axis 
of -Bft ec is uncorrelated with the thrust axis of -B^ag- In a qq event, the decay products lie 
in the two jets which are back-to-back, so the two thrust axes tend to be collinear. Before 
training, a loose requirement of | cos^tb^oI < 0.9 is imposed that retains 90% of the signal 
while rejecting 50% of the continuum background. The Fisher discriminant is also required 
to satisfy —3 < J^bl/qq < 2. 

Background from charm (b — > c) decays are found to be negligible and are not considered, 
while charmless (b — > u, d, s) decays of the B meson do contribute to candidates in the 
kinematically studied region, though rarely to the signal region. 

On average, 1.01 B candidates are reconstructed per signal event and the fraction of 
mis-reconstructed events is negligible. On the rare occasion that a single event contains 
multiple candidates, the B candidate having the two track with the highest momentum in 
the lab frame is selected. 

Since the -Br cc and B^ ag mesons are approximately at rest in the T(4S) CMS, the differ- 
ence in decay time between the two B mesons, At, can be approximately determined from 
the displacement in z between the final state decay vertices, 

At ~ (^ec ~ ^Tag) = 

m ~ /3 7 c ~ he { ) 

The vertex of reconstructed B candidates is determined from the charged daughters, 
with a further constraint to lie within the known IP. The IP profile is smeared in the plane 
perpendicular to the z-axis to account for the finite flight length of the B meson in this 
plane. To obtain the At distribution, we reconstruct the tag side vertex from the tracks not 



used to reconstruct B^ ec [19J. The events must satisfy the requirements |At| < 70 ps and 
^-Rec.Tag < 500, where /iR ec ,Tag is the multi-track vertex goodness-of-fit, calculated in three- 
dimensional space without using the interaction-region profile's constraint [3|. To reduce 
the necessity of also modelling the event- dependent observables describing the variable At 
resolution in the fit, the vertex uncertainty is required to be <j^ cc,Tag < 200 /im for multi-track 
vertices and cr^ ec,Tag < 500 /^m for single-track vertices. 

The flavor tagging routine is employed as described in Ref. [2l[ . The tagging information 
is represented by two parameters, the B^ ag flavor q and the tagging quality r. The parameter 
r is an event-by-event MC determined flavor tagging quality parameter that ranges from 
r = for no flavor discrimination to r = 1 for unambiguous flavor assignment. Due to 
the non-zero probability of mistagging w, the CP asymmetry is diluted by a factor 1 — 
2w. The measure of the performance of the flavor tagging algorithm is the total effective 
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tagging efficiency e c g = 6Ta g (l — 2w) 2 , as the statistical significance of the CP parameters is 
proportional to (1 — 2to)^/eTag, where ex a g is the raw tagging efficiency. These are determined 
to be e eff = 0.284 ± 0.010 and e cS = 0.301 ± 0.004 for SVD1 and SVD2 data, respectively. 

Differences from the previous Belle analysis include a new tracking algorithm and the 
inclusion of the event shape Fbh/qqi the fit rather than optimizing selection criteria for 
this variable. As this results in a large increase of the continuum background, a reduced 
kinematic fit region in M bc and AE is applied in order to reduce this background without 
significant loss of signal events. These changes increase the detection efficiency by 19% over 
the previous analysis according to MC simulation. 

IV. EVENT MODEL 

The CP violation parameters are extracted from a seven- dimensional extended unbinned 
maximum likelihood fit to M bc , AE, T h ij q q, At and q from a data sample divided into 

7 bins (/ = 0..6) in the flavor-tag quality r and 2 SVD configurations s. Seven categories 
are considered in the event model: B° — > -k + ti^ signal, B° — > K + n~ , B° — > K~n + and 
B° — > K + K~ peaking backgrounds, continuum, charmless neutral and charmless charged 
BB decays. In most categories, the linear correlations between fit variables are small, so the 
probability density function (PDF) for each category j, is taken as the product of individual 
PDFs for each variable P-' s (M bc , AE, F bl/qq -, £+ /?r , C~ /n , At, q) = V l ' s {M bc ) x V l ' s (AE) x 
^P l ' s {-Fbb/qq) x ^ :S (^k/7t' ^k/tx) x P l,s {At,q) in each l,s bin, unless otherwise mentioned. 

A. Peaking Models 

The four peaking shapes including the signal are determined from reconstructed MC 
events. The PDFs for M bc and AE are taken to be the sum of three Gaussian functions 
where the tail Gaussians are parameterized relative to the core, incorporating calibration 
factors that correct for the difference between data and MC. These factors calibrate the 
mean and width of the core Gaussian and are determined from a large-statistics control 
sample B + — > D°[K + h~]it + . The PDF for T^ij^ here and also for the two background 
BB categories are the sum of three Gaussians in each flavor-tag bin I. The shape for the 
four peaking channels are identical and the BB background shapes are fixed from the signal 
model except for the mean and width of the core distribution. Calibration factors that 
correct for the shape difference between data and MC are similarly incorporated into the 
core mean and width. The shape is modelled with a two-dimensional histogram that 
has been corrected for the difference between data and MC in PID as determined from an 
independent internal study with inclusive D* + — > D°[K~ 7r + }n^ ow decays. The At,q PDF for 
B° — > 7r + 7r~is given by 

V l ^{At,q) = 4rgQ |l-gA^ + g(l-2^ s )x 

Acp cos Am^At + Sep sin Am^At 

which accounts for CP dilution from the probability of incorrect flavor tagging w l,s and 
the wrong tag difference Aw l,s between B° and B°, both of which are determined from 



j ®R S B0B0 (At), (4) 
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The physics parameters rgo and Am^ are fixed to their 
This PDF is convolved with the At resolution function 



for neutral B particles R S B0S0 , described in Ref. [19(. We consider the At,q distributions for 
the flavor-specific B° — > K + ir~ and B° — > K~~ir + peaking backgrounds separately with 



-\At\/T g0 r \ 

V l ^{At,q) = — 1 1 - qAw l > s T g(l - 2w l ' s ) cos Am d At j R% 0S0 (At). (5) 

For the B° — >• K + K~ peaking background, the PDF for At,q is the same as for B° — > 7r + 7r~ 
signal, however, the CP parameters are fixed to zero. To account for the broad underlying 
At events not described by either peaking or background PDFs, a broad Gaussian outlier 
PDF is introduced for every category 

V l 6 s ut (At,q) = ±G(At;0,a s o J. (6) 



B. Continuum Model 

The parameterization of the continuum model is chosen based on data taken 60 MeV 
below the T(45) resonance; however, all the shape parameters of Mb c , AE, T b i/ q q and J~-^/ n 
are free in the fit to extract the CP parameters. Since continuum is by far the dominant 
component, extra care must be taken to ensure this background shape is understood, so 
correlations above 2% are considered. The PDF for M^ c is an empirical ARGUS function [23], 
while AE is parameterized by a straight line in each flavor-tag bin with a slope depending 
linearly on jF^/ qq 

P%(AE\JF a/99 ) = 1 + (ri s + p{T bl/qq )AE. (7) 
The Fbl/qq shape is observed to shift depending on PID region, so the PDF is a sum of two 
Gaussian functions in two PID regions, C^, n <= 0.5 and (£>K/ n or ^K/n) > ^ sman cor_ 
relation between the C K ,^ shape and flavor-tag q is also observed due to the ss component of 
continuum. For example, if a B^ cc candidate is reconstructed with a single high-momentum 
K + , it inhabits the flavor specific K + ir~ sector of The accompanying K~ can then 

be used as part of the flavor tagging routine which will lead to a preferred flavor-tag of B°. 
Thus, the distribution is enhanced in the K + n~ region and is depleted in the K~n + 

region for q — — 1. To account for this effect, we model with an effective asymmetry 

A l q y q, that modifies the full two-dimensional PID histogram model H l ' s {C^-,^ C^-,^), in each 
/, s bin depending on the flavor-tag, 



where 



= -«o1^ /7r - C+fj* if C K/7r - C + K/7T < 0, (9) 
which we hereafter refer to as the "manta ray" function. The At model 



-\At\/r q9 

h)— + fs6(At-^) 



RiJAt), (10) 
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contains a lifetime and prompt component to account for the charm and charmless contri- 
butions, respectively. It is convolved with a sum of two Gaussians 

R S qq(^t) = (1 — /tail)^(^ Mmeam ^main ") + fLi\G(At] /Xycain ^main^tail) > (H) 

which uses the event-dependent At error constructed from the estimated vertex resolution 
cr = ( <J o"^ cc + °"Tag) / /^7 C as a sca l e factor to the width parameters S^ain aric ^ ^taii- 

C. Model 

The charmless BB background shape is determined from a large sample of MC containing 
generic b — > u,d,s transitions and is subdivided into neutral and charged B samples. A 
sizeable correlation is seen between M^ c and AE and is taken into account with a two- 
dimensional histogram. A similar correlation between the flavor-tag and , similar to 

that in continuum is also seen here. However, due to B°B° mixing in the neutral BB 
background, this effect must be modelled along with At, q. For the neutral BB background, 
the PDF is given by 

^o 50 (£^,At,g) = tf^^/Jx 

e -\M\/T B0B0 



4r 



BOB 



jl + qA% S0 (C+ /w ,C^ /w ) cosAm d Atj <g> R% 0B0 (At), 

(12) 



and the charged BB background PDF is given by 

T> l > s (C± At n\- 1 + qA B+B-( C K/m C K/J rrl,s (r + n- eT^^B- 

^B+B-y L K/^ At ^J ~ 2 y L K/^ L K/J 2t B+B - ® R B+B-{ At )i 

(13) 

where A B S g are manta ray functions for each BB category and Rb+b- is the At resolution 
function for charged B particles. Since reconstructed BB events may borrow a track from 
the tag side, the average At lifetime tends to be smaller and is taken into account with the 
effective lifetime, t bb . 



D. Full Model 

The total likelihood for 559797 B° — > h + h~ candidates in the fit region is 

■ £j N f 12i, B fj' s Nl ' s 



e 



r v .v;/;-^;-(.\ca/-. ^ d* c-* At\ q % (i4) 



]\f | 11 i J j r 3 \ be i bb/qq 1 K/'K 1 K/m 

l,s i=l j 

which iterates over % events, j categories, I flavor-tag bins and s detector configurations. The 
fraction of events in each flavor-tag bin I, for category j, is denoted by /y s . The fraction 

of signal events in each l,s bin, / g '/, has been calibrated with the B + — > D [K +/ k~]7t + 
control sample. Free parameters include the B° — > K + n~ and B° — > K + K~ yields, iVf= and 
N s BOm \ the remaining N s BOm yields are fixed to N*£%l = (0.269±0.010)A^yj > 1 and = 
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(0.268 + 0. 004) A^f]^ 2 as determined from MC. In addition, the signal M^ c calibration factors 
that account for the difference between data and MC are free in the fit to data as well as 
all shape parameters of the continuum model with the exception of the At parameters. In 
total, there are 116 free parameters in the fit. 

In contrast to the previous Belle analysis (BJ, here all variables are fit simultaneously to 
determine the component yields and CP violation parameters. The previous analysis applies 
a 2-step procedure where the event-dependent component probabilities are calculated from 
a fit without At and q. These are then used as input in a fit to At and q to set the fractions 
of each component to determine the CP parameters. Our procedure has the added benefit 
of further discrimination against continuum with the At variable and makes the treatment 
of systematic uncertainties more straightforward at a cost of CPU and analysis complexity. 



V. FIT RESULT 

A fit is performed to the data, with the projections shown in Figs. |3] and HI and the 
following CP violation parameters are obtained 

A C p(B° -»■ 7T+7T-) = +0.33 ± 0.06 (stat) ± 0.03 (syst), 

S C p(B° -»■ 7T+7T-) = -0.64 ± 0.08 (stat) ± 0.03 (syst), (15) 

where the first uncertainty is statistical and the second is the systematic error, which will be 
discussed in the next section. Our results confirm CP violation in this channel as reported 
in previous measurements and other experiments j5|-0| , and they are in marginal agreement 
with the previous Belle measurement at the 1.8<r level. Besides the increase in data sample, 
the difference between results can also be attributed to the new tracking algorithm and the 
19% increase in detection efficiency These are also the world's most precise measurements 
of time-dependent CP violation parameters in B° — > 7t + tt~. The statistical correlation 
coefficients between the obtained parameters is +0.10. The peaking event yields including 
signal are N(B° ->■ tt+tt-) = 2964 ± 88, N(B° -> K+tt') = 9205 ± 124 and N(B° -> 
K + K~) = 23 + 35 where the uncertainties are statistical only. From the obtained yields in 
the fit, the relative contributions of each component are found to be 0.5% for B° — > tt + tt^, 
1.7% for B° —7- A' + 7r~, 97.6% for continuum and 0.2% for BB background. 

Using Eq. [2] and input from other Belle publications (H], 25|| , an isospin analysis is per- 



formed to constrain <p2- A x 2 is constructed for the 5 amplitudes shown in Fig. [2] accounting 
for the correlations between our measured physics observables used as input. The \ 2 is then 
converted into a confidence level (CL) as shown in Fig. [5j The region 23.8° < <p2 < 66.8°, 
is ruled out and the constraint on the shift in <p 2 caused by the penguin contribution is 
|A0 2 | < 44.8° at the la level including systematic uncertainties. 



VI. SYSTEMATIC UNCERTAINTIES 

Systematic errors from various sources are considered and estimated with independent 
internal studies and cross-checks. These are summarized in Table [TT1 Uncertainties affecting 
the vertex reconstruction include the IP profile, charged track selection based on track helix 
errors, helix parameter corrections, At and vertex goodness-of-fit selection, Az bias and SVD 
misalignment. The fit model uncertainties including the fixed physics parameters r^o and 
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FIG. 3: Projections of the fit to the data enhanced in the B° — > tt + it signal region. The points 
with error bars represent the data and the solid black curves or histograms represent the fit results. 
The enhancements, M bc > 5.27 GeV/c 2 , \AE\ < 0.04 GeV, T h i jm > 0, £± /w < 0.4 and r > 0.5, 
except for the enhancement of the dimension being plotted are applied to each projection, (a) (b) 
(c) (d) and (e) show the Mb c , AE, C^,, and F^b/qq projections, respectively. Blue hatched 

curves show the B° — > tt + it~ signal component while green dotted curves show the B° —> K^-k^ 
peaking background component. Dashed red curves indicate the total background while the purple 
dash-dotted curve shows the BB background component. 
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FIG. 4: Background subtracted time-dependent fit results for B° — > tt + it~. (a) shows the At 
distribution for each -B^ ag flavor q. The solid blue and dashed red curves represent the At dis- 
tributions for B° and B° tags, respectively, (b) shows the asymmetry of the plot above them, 
(N B o — N^o)/(N B o + Ngo), where N B o {N B o) is the measured signal yield of B° (B°) events in 
each bin of At. 




FIG. 5: Difference 1-CL, plotted for a range of 4>2 (a) and (A^l (b) values as shown by the solid 
curve. The dashed line indicate the la exclusion level. 

Arrid, parameters describing the difference between data and MC, At resolution function 
parameters, as well as the flavor tagging performance parameters w and Aw, are varied 
by ±lcr. The parametric and non-parametric shapes describing the background are varied 
within their uncertainties. For non-parameteric shapes (i.e. histograms) we vary the contents 
of the histogram bins within ±lo\ The fit bias is determined from pseudo-experiments by 
searching for a difference between the generated and fitted physics parameters. Finally, a 
large number of MC pseudo-experiments are generated and an ensemble test is performed 
to obtain possible systematic biases from interference on the tag-side arising between the 
CKM favored bd —> (cud)d and doubly-CKM-suppressed bd —> {ucd)d amplitudes in final 
states used for flavor tagging j26[. 
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TABLE II: Systematic uncertainties of the measured physics parameters. 



Category 


SAcp(tt + tt-) (10- 2 ) 


SScp(tt + tt-) (10- 2 ) 


IP profile 


0.1 


1.2 


track selection 


0.3 


0.3 


Track helix errors 


0.0 


0.0 


At selection 


0.0 


0.0 


Vertex quality selection 


0.4 


0.2 


Az bias 


0.5 


0.4 


Misalignment 


0.4 


0.2 


Tgo and Am<j 


0.1 


0.1 


Data/MC shape 


0.2 


0.2 


At resolution function 


0.8 


2.0 


Flavor tagging 


0.4 


0.3 


Parametric shape 


0.2 


0.3 


Non-parametric shape 


0.4 


0.6 


Fit bias 


0.6 


0.9 


Tag-side interference 


3.2 


0.2 


Total 


3.5 


2.7 



VII. CONCLUSION 

We present an improved measurement of the CP violation parameters of B° — > -k + ti^ 
decays, 

Acp(B° tt+tt - ) = +0.33 ± 0.06 (stat) ± 0.03 (syst), 

S C p(B° ->■ 7r + 7r~) = -0.64 ± 0.08 (stat) ± 0.03 (syst), (16) 

confirming CP violation in these channels reported in previous measurements and other 
experiments jH|-0]. These results from the full Belle data sample after reprocessing with a 
new tracking algorithm and with an optimized analysis performed with a single simultaneous 
fit, supersede those of the previous Belle analysis |5|. They are now the world's most precise 
measurement of time-dependent CP violation parameters in B° — > 7t + tt~, ruling out the 
range 23.8° < <p 2 < 66.8°, at the la level. 
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